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Abstract. Using the conformally invariant Cotton tensor, we define a geometric flow, 
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initial metrics into conformally flat ones, and is somewhat orthogonal to the Yamabc 
flow, the latter being a flow within a conformal class. We define an entropy functional, 
and study the flow of nine homogeneous spaces both numerically and analytically. In 
particular, we show that the arbitrarily deformed homogeneous 3-sphere flows into the 
round 3-sphere. Two of the nine homogeneous geometries, which are degenerated by 
the Ricci flow, are left intact by the Cotton flow. 
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1. Introduction 

In this paper, inspired by Hamilton's work on the Ricci flow [lj that culminated in 
Perelman's proof of Thurston's geometrization conjecture and the longstanding Poincare 
conjecture (2j [3j, we introduce a new flow, based on the conformally invariant Cotton 
tensor. This flow is exclusive to three dimensions and its fixed points are the conformally 
flat metrics. 

Historically, the Ricci flow equations seem to have first arisen in a work by Friedan 
[I] that deals with the renormalization group (RG) flow in a 2-dimensional nonlinear 
sigma model. Just like any other coupling in quantum field theories, the Riemannian 
metric, considered as a coupling 'constant' between the sigma model fields defined on a 
manifold M, runs as the energy scale changes: 

(3 g = d t gij = ka'Rij, (1) 

where k is a constant uniquely fixed by the 1-loop result, and a' is a dimensionful 
parameter (such as the tension). Of course, in general, the beta function for the RG 
flow of the metric will have infinitely many terms, including the Ricci tensor. In fact, in 
[I], the 2- loop result with a Riemann-square correction is given. In string theory, which 
is again a two dimensional sigma model with additional fields, gravitation (including 
General Relativity as the lowest order approximation) that governs the dynamics of 
the target space, is derived by the requirement that the beta function vanishes as a 
result of Weyl invar iance on the worldsheet. Since a non-anomalous Weyl symmetry 
is crucial for string theory, in the physics literature, it is common practice to set 
P g = and consequently, end up not with a flow, but with a usual gravity theory 
(containing perhaps some additional fields and higher curvature terms). On the other 
hand, Hamilton pQ introduced the Ricci flow as a means to deform a given metric on 
a manifold and developed a programme aiming to prove the Poincare conjecture. In 
particular, he classified closed 3-manifolds of positive Ricci curvature. 

Obviously, one could introduce many different flows, or depending on one's point 
of view, derive flows that come from the beta functions of two dimensional non-linear 
sigma models. As examples of the former, see (5j El [7j. For an example of the latter, see 
[8] where a string inspired flow with scalar, Maxwell and anti-symmetric tensor fields is 
studied. 

Our motivation was to find a flow whose critical points are the conformally flat 
spaces, and this is achieved in three dimensions with a unique 2-tensor called the Cotton 
tensor. We study the behaviour of the nine homogeneous geometries under this flow in 
detail. Four of these are fixed points of our flow. [We remark that the eight geometries 
appearing in Thurston's geometrization conjecture are maximally symmetric; viewed in 
this way, two of the nine (I som(M. 2 ) and R 3 ) collapse to the same maximally symmetric 
geometry] 

Interesting physics could come out of the Ricci flow or other flows such as the one 
we introduce here, especially in Euclidean quantum gravity defined by the path integral. 
This is somewhat an unexplored territory. Yet, we can point out several works in this 
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direction: In [5], the Ricci flow is applied to four dimensional black hole physics as a 
means to understand the phase structure of this theory in Euclidean space. The Ricci 
flow is also used [TDJ to find inequalities in General Relativity regarding the evolution 
of the area of a surface and the enclosed Hawking mass. See [11] for a compilation of 
possible applications of the Ricci flow in physics. 

The outline of the paper is as follows: In section 2, we introduce the Cotton flow 
and show how some geometric quantities evolve accordingly. In section 3, we give an 
entropy functional which is non- decreasing under the Cotton flow. This also gives us a 
gradient formulation of the Cotton flow. In section 4, we show that the Cotton flow is in 
a certain sense orthogonal to the Yamabe flow. Section 5 is the bulk of our paper where 
we study the evolution of homogeneous spaces under the Cotton flow both numerically 
and analytically. In the conclusion, we summarize our results and mention some topics 
for further study. Details of the computations regarding the homogeneous geometries 
are given in the appendices. 



2. The Cotton flow 



In dimensions greater than or equal to 4, the Weyl tensor determines whether or 
not a given manifold is locally conformally flat. In 3-dimensions, the Weyl tensor 
vanishes identically, and the role of the Weyl tensor is played by the Cotton tensor 

H21 MM El EE! 

Amn i 

C7a = —V m (# B (2) 

where e is a tensor density, in an orthonormal frame e 123 = 1, Rij is the Ricci tensor, 
and 1Z the curvature scalar. [Note that the Cotton tensor is also referred to as the 
'Cotton- York tensor' since York reintroduced [15] and used it extensively in the initial 
value formulation of General Relativity] Here and throughout, the signature of the 
space will be (+, +, +), the indices will range from {1, 2, 3}, and the sign convention for 
the Riemann tensor will follow from [V,, Vjjffc = Rijk i>j. The tensor C u is covariantly 
conserved, symmetric and traceless (see [17], where the Cotton tensor was used to 
introduce the "topologically massive gravity"). The symmetry can be directly seen 
from the following representation of C JJ : 

= U mn V m R j n + e^VmK) . (3) 

The tensor C % j is conformally invariant with weight 0. 

Let M be a smooth 3- manifold with a positive definite metric gij. We will consider 
the flow defined by the equation 

d t gij = nCij, (4) 

where k is a positive constant (note that t is not one of the coordinates of the manifold). 
Without loss of generality, we may scale t to set k = 1. The reason for taking k > is 
to make homogeneous metrics on the 3-sphere converge to the round sphere rather than 
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diverge from it. Since the right hand side is traceless, yjg is independent of t although 
it could depend on the coordinates. In contrast, the Ricci tensor is not traceless in 
general, so a normalization factor is necessary for the Ricci flow (OQ) in order for the 
volume density or volume to be preserved. 

Assuming that g^ flows under the evolution equation (jlj), one can compute how 
some other geometric quantities evolve: 

d t r Jk = \g il {VjC kl + v k Cfl - v t c Jk ), 

dtRij = 3R[(iCj) 1 — R lm Ci m gij — -TZCij — -V 2 CV,-, 

d t n = -c^Rij, 

dtC ij = e mk{i V ^ dtRj ) k + e mk{i dtVj ) ^ fe) 

d t {R, k w k ) = m 3l c l k w k - mR jk k - R 3k v 2 c, k , (5) 

where the usual shorthand notation for symmetrization is used: A(i\j\ k ) = (Aij k +A k ji)/2. 

The Cotton flow equations form a system of third order, nonlinear partial 
differential equations. About a flat background, the linearization of the system, setting 
gij = o~ij + hij, gives 

t h ij = ~^d 2 d k {e M hj + e jkl hf) + + (6) 

where 

C = -\e tkl d k d n h nl . (7) 

and d 2 denotes the Laplacian in flat space. The £ terms can be viewed as coming 
from a diffeomorphism of the manifold. Unlike the Ricci flow case, whose linearized 
form is of heat equation type in the highest order after DeTurck's modification by a 
diffeomorphism |18j . our equation is of third order. In this respect, it seems impossible to 
apply results from elliptic operator theory directly for the short time existence problem. 
Understanding the conditions for the initial metric under which the equations are locally 
well-posed seems to be a delicate problem, which will not be studied in this paper. 

3. Cotton entropy 

In order to show that the Ricci flow can be regarded as a gradient flow, Perelman [2] 
defined an entropy functional T{g, f) = j M (R + \ V f\ 2 )e~f ^/gd 3 x, where g denotes the 
metric, and / is a scalar field. He proved that if d t f = d t In ^fg and g obeys the Ricci 
flow equations, then JF is a (non-strictly) increasing functional. For the Cotton flow, an 
analogous functional is the well known gravitational Chern-Simons action |17j : 

Hg) = -\j M t l3k T^m(d,T m kl + h m jn T n ki y 3 x. (8) 

Note the absence of an extra scalar field, in contrast to the case of the Ricci flow. 
The functional T is conformally invariant, and up to a boundary term, diffeomorphism 
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invariant. Defining d t gij = fy, after a straightforward but lengthy computation, one 
obtains 

IF = / M vW'^.r. ( 9 ) 

The choice i>y = CV,- gives the steepest descent, and leads to the Cotton flow d t gij = Cij. 
The functional T{g) is increasing if g evolves under the Cotton flow. The Cotton entropy 
is constant if and only if Cy is identically zero, which means that M is locally 
conformally flat at all points. 

4. Cotton flow and the Yamabe flow 

One can define various flows in the space of metrics. One of these is the Yamabe flow 
defined by the equation [6J 

9 T gij = -Tlgij- (10) 

Under this flow, the conformal class of a metric does not change. The Yamabe flow 
was introduced in order to solve Yamabe's conjecture [Bj, which states that any metric 
is conformally equivalent to a metric with constant scalar curvature. The Cotton flow 
has a somewhat complementary behaviour since by definition it changes the conformal 
class unless the metric is conformally flat. We will now make this behaviour precise in 
a sense, by computing the commutator of the vector fields for these two flows. 

Let t and r denote, respectively, the parameters for the Cotton and the Yamabe 
flows. Then 

[d t , d T ]gij = d t d T g i: j - d T d t g i:j = -d t (Kgij) - d T C i:j . (11) 

Since C % j is conformally invariant, d T C % j = 0. From here, it follows that <9 r CV, = —IZCij. 
Using the third equality of (jSJ), one gets d t {JZgij) = IZCij — R kl Ckigij. Therefore, 

[d t) d T }g l3 = R hl C kl g ir (12) 

Since the result is proportional to the metric, the commutator of the two flows gives 
another flow preserving the conformal class, just like the Yamabe flow. 

5. Cotton flow on homogeneous 3-manifolds 

In this section, we will be interested in studying the behaviour of the Cotton flow on 
homogeneous three manifolds. There are nine such homogeneous geometries: 577(2), 
SL(2, R), IsomiM 2 ), Solv, Nil, M 3 , H 3 , S 2 xR, H 2 xR. As explained in the introduction, 
they give rise to the eight maximally symmetric geometries which are the basic building 
blocks that take the stage in Thurston's geometrization conjecture for three manifolds 
[T9] . This conjecture, which includes the Poincare conjecture as a special case, was 
settled by Perelman's work (21 [3] . 

The first six of these nine geometries can be obtained by the following construction: 
Suppose that G is a 3-dimensional, unimodular, simply connected Lie group, and that 
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g is a left-invariant metric on G. Let M be a compact 3-manifold obtained as the 
quotient of G with respect to a discrete subgroup. The analysis of the Ricci flow on 
these manifolds was carried out in detail in [201 EJJ (see also [22] ) . [A similar analysis 
for cross curvature flow was carried out in [23].] As in the case of the Ricci flow, the 
Cotton flow equations on a homogeneous 3-manifold reduce to a set of three coupled 
ordinary, nonlinear, autonomous differential equations. 

We now proceed to find this system of differential equations. In [23], it is 
shown that there exists an orthogonal left-invariant frame {F±, F2, F3} on G such that 
[F X ,F 2 ] = 2vF 3 , [F 2 ,F 3 ] = 2AFi, [F 3 ,Fx\ = 2/iF 2 where /z, v, X G {-1,0,1} (see also 
|22j). Let {a; 1 , uj 2 , uj 3 } be the dual basis of 1-forms to {F\, F2, F3}. In this basis, the 
metric g takes the form 

g = Aitp 1 ® uj 1 + B(t)uj 2 ® uj 2 + C{t)uj 3 <g> u; 3 . (13) 

Defining e 1 = v^Lu 1 , e 2 = \I~Buj 2 , e 3 = \fCuj 3 we obtain an orthonormal coframe, with 
respect to which the metric takes the form g = e 1 <S> e 1 + e 2 <S> e 2 + e 3 <8> e 3 . One has 

gL/ = 2Ac<j 2 A u 3 , du 2 = 2fiu 3 Aw 1 , dw 3 = 2w ! A w 2 . (14) 

The Cotton flow equations can be written as dte a = *C a , where C a is the Cotton 2-form 
and * denotes the Hodge dual. Since the Cotton flow preserves the volume density, from 
now on we assume that ABC = 1. Referring the reader to the appendix for the details 
of the computation, we state the equations which take the form 

A A 

— = 4A[-X 2 A 2 (fiB + uC- 2XA) - (fiB + uC)(fiB - vC)\ 

^ = AB[-fi 2 B 2 (uC + XA- 2fiB) - [yC + XA){yC - XA) 2 }, 
dC 

= AC[-u 2 C 2 (XA + 12B- 2vC) - (XA + /iB)(XA - fiB) 2 }. (15) 

(JjL 

The entropy functional T defined in section |3] can be alternatively computed using 

which for our case yields 

jr = y [AXfxu-2(XA + fiB-uC)(XA-fxB + uC)(-XA+fiB + uC)]*l,(16) 

where *1 = e 1 A e 2 A e 3 . Recall that by the result in section [31 T is a non- 
decreasing functional. [The name "entropy" may be misleading since T is not necessarily 
nonnegative. In fact, the maximum value of T in the SU(2) case below is —2. All 
that matters is that T is non-decreasing.] The following simpler functional, which was 
obtained by manipulating ( Tl5i) . also acts like an entropy in certain cases, and will also 
be useful in what follows: 

1 /JI 2 I/ 2 \2„2 \2.,2 N 
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Note also that, using (TT5|) 



2 2 \2 2 

' ^ V (fj,B + vC)(nB - uC) 2 + ^-(vC + XA)(uC - XAf 



dt A 2 ^ ' B 2 

C 



\ 2 2 

+ -^(\A + fiB)(\A-nB) 2 . 



5.1. SU(2) 



When we take A = \i — v — — 1, the Lie group G is isomorphic to SU(2), which is 
homeomorphic to S 3 as a topological manifold. The proof that the flow equations (|T5l) 
take any initial homogeneous metric to the round metric A = B = C = 1 on S 3 is 
remarkably simple. Since A = fi = v = —1, one has f(t) > 3/8, df /dt < 0, and in both 
cases equality occurs if and only if A — B — C — 1. This proves the claim. 

It seems impossible to find an exact analytical solution to this system. Therefore 
we carry out an estimate analysis similar to that in [20]. From (I15p . one gets 



d<yA B "> = 8(B - A)[2(B 3 + B 2 A + BA 2 + A 3 ) - C(A + B) 2 - C 3 ].(19) 

(JjL 

Since the right hand side vanishes when A = B, by the uniqueness of solutions, the 
ordering between A and B, and likewise the ordering of the three variables cannot 
change during the flow. Therefore, without loss of generality, we may assume that 
A > B > C due to the symmetry of the equations. From the last equation of (1151) . it is 
clear that C is not decreasing, so A > B > C > Co where Cq is the initial value of C. 
Moreover, 

d ^ A ~ °^ = -8{A - C)[2C 3 + 2C 2 A + 2CA 2 + 2A 3 - B{A + C) 2 - B 3 ] 

(aiL 

< -8(A - C)[2C 3 + 2C 2 A + 2CA 2 + 2A 3 - A(A + C) 2 - A 3 } 

< -2AC 3 (A-C). (20) 

We deduce that A — C < (A — C ) exp(— 24Cgt). Since the difference between the 
greatest and smallest of A, B, C decays at least at this rate to 0, (A, B, C) also converges 
at least at this rate to (1, 1, 1). 

Linearizing the system around the fixed point (A, B) = (1,1) after setting C = 
1/AB gives 

dA „ dB , , 

M = ~ 24A - M-- 2iB - < 21) 
which is in accordance with the estimate above since max(Co) = 1. 

One can alternatively solve the equations numerically. Starting from various initial 

values for (A , B , Co), a sketch of the solution curves can be seen in figure 1. 

5.2. SL(2,R) 

When we take A = /i = — 1, u = 1, the group G is isomorphic to the universal cover of 
SX(2,R). From ([15]) one gets 

d<yA ~ B "> = 4{B-A)[2{B 3 + B 2 A + BA 2 + A 3 ) + C{B 2 + A 2 )+2 + C 3 }.{22) 

(J/L 
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Figure 1. The solid lines indicate flows, for various initial metrics on S 3 : (Aq,Bo) = 
(0.2,0.8), (0.5,0.6), (0.7,0.3), (1,1.9), (1.6,0.9). The horizontal axis is for A(t) and 
the vertical axis is for B(t). C(t) is not depicted as it runs according to C = 1/AB. 



Since the right hand side is when A = B, we deduce that the ordering between A 
and B does not change throughout the flow. Without loss of generality, assume A > B 
(note that the equations for A and B are symmetric). The equation for A reads 

A A 

= 4A[A 2 B - A 2 C - 2A 3 + B 3 + B 2 C - BC 2 - C 3 ] 

< —A(ABC 2 + AC 3 ) = -4(C + AC 3 ). (23) 

This shows that A is strictly decreasing (However B is not necessarily decreasing). The 
equation for C gives a lower bound 
dC 

^- = AC[C 2 (A + B + 2C) + {A + B)(A- B) 2 ] > 8C\ (24) 

Therefore, C — > oo in a finite time tf. We want to show that A(tf) = B(tf) = 0. Since 
there exists a point ti of time after which C > A + B, we have dC/dt < 16C 4 . For 
t > k 

< _4^c<4 < o (25 ) 
dt ~ ~ ( J 

which implies that A(tf) = 0. Since B{tf) < A(tf), we also get B(tj) = 0. In the 

terminology of [20J, the geometry always tends to a "cigar degeneracy" in finite time. 

Some numerical solutions of the system are given in figure 2. 

5.3. Isom(R 2 ) 

When we take X — fi — —1, u — 0, G is isomorphic to the universal cover of the group 
of Euclidean motions of IR 2 . From (TT5T) one obtains 



d ^ A B) =8(B 4 -A 4 ), ^=4C(A 2 -B 2 )(A-B). (26) 
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Figure 2. The solid lines indicate flows, for various initial metrics on a homogeneous 
geometry modelled on SX(2,R): (A ,B ) = (0.5,1), (2,0.7), (1.5,1.8), (0.2,2). 



From these two equations, one gets 

d(A-B 



C(A-B) 



+ 2[(B-A) 2 + l)^ = 0. 



(27) 



dt ' ~V ' G) dt 

The function C 3 is an integrating factor for this equation. So, there is a second conserved 
quantity besides ABC = 1: 



C\A - B) 



k. 



(28) 



where k is a constant. Since there are two conserved quantities and three variables, 
the phase space foliates into 1-dimensional leaves. From (1261) . we see that the ordering 
between the symmetric variables A and B doesn't change during the flow. Assume 
without loss of generality that A > B. Then, 



d A 

^ = 4A(A 2 B - 2A 3 + B 3 ) < 0, 

dB 
~dt 



AB(A 3 + AB 2 - 2B 3 ) > 0. 



(29) 
(30) 



Thus A is decreasing, B is increasing, and they are both bounded away from 0. This 
implies that the flow exists for all t. If Bo is the initial value of B, then since B > Bq, 
from (1261) one gets 

d{A - B) 



dt 



= -8(A - B)(B 3 + B 2 A + BA 2 + A 3 ) 
< -32B 3 (A - B). 



(31) 
(32) 



So, A — B decays to 0. In fact, using the conserved quantities, the equations can be 
integrated. From ([28]), it follows that (A - B) 2 = ^ - ^. Using AB = 1/C, one 
obtains (A + B) 2 - k 

dC 



3C 

+ Then, the second equation in (12"6"1) yields 



dt 



4( 



C 3 



k 4 
C 1+ 3C 
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Defining 



one finds 



C 



(?(«'-*) 



1/3 



2 — 



32 



3A; - 2m 2 dt 

whose integration leads to an implicit definition of u(t) as 



Vk 



u 

~2 + 2V6 



tanh 



32 



and I is another constant. The remaining metric functions follow easily as 

-2/3 

k)' 



A 



1/3 2 

2U (W 
1 /3 



+ V3k - 2u 2 



-2/3 



(m 2 - k) J (« - V3A; - 2u 2 ) 



5 , . 
2 V4 

Since A - B — ► 0, by (|2gjl . the limiting value of C as t -> oo is (3fc/8) 1/3 . In the limit, 
A = B and the geometry becomes flat. Even though we have analytically solved the 
equations, we also present some numerical solutions of the system for various initial 
conditions in figure 3. 



Metric Components 




Figure 3. The solid lines indicate flows, for various initial metrics on a homogenous 
geometry modelled on Isom(R 2 ): {A ,B ) = (2,0.2), (1,2), (1.6,0.5), (0.7,1.8), 
(0.4,1). 



The conserved quantity (1281) can be, more invariantly, written in terms of the Ricci 
tensor and the scalar curvature as 
16 4 TZ 3 Rtc 2 

const. (33) 



6 {?>1Z 2 -Ric 2 Y 

To get a unitless quantity, one should divide this by Vol(M) 2 , if the manifold is of finite 
volume. 
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When we take A = — 1, \i — 1, v = 0, the group G is isomorphic to the group of 
isometries of the Minkowski plane. The equations ( fT5i) take the form: 
dA . , , .s dB 

~dt 

dC 



AA{A 2 {B + 2A) + B 3 



dt 



AB(B 2 (A + 2B) + A 3 ), 



— = -AC(B 2 - A 2 )(A + B). 
dt 



(34) 



The first two equations imply that A is decreasing and B is increasing. In fact, 
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Figure 4. The solid lines indicate flows, for various initial metrics on a Solv geometry: 
(A ,B ) = (2,0.1), (1.5,0.3), (1.8,1). 



dB/dt > 8B 4 , which shows that B — > oo in finite time. As in section 15.31 one can 
find a second conserved quantity which now reads 

(A + B) 2 C A - -C 3 = k, (35) 
3 

where k is a constant. Then, since A is decreasing, C must decay to at a rate like 
1/\^B. Recalling ABC = 1, we see that A also decays to like 1/VB. All geometries 
in this class thus tend to a cigar degeneracy. Again, the equations can be completely 
integrated to give 

1/3 \ -2/3 , 

A = -(-(k-u 2 )) (V3k - 2u 2 - u) , 



where 



2 V4 

1/3 \ -2/3 , /3 \ 1/3 

B = ^(k-u 2 )) (V3k—2^ + u), C=(^k-u 2 )) 



u Vk , i . j 7— . 32 



tanh _1 (uV2/3fc) = — t + I 



2 2VQ 3 
and I is another constant. We present some numerical solutions of the system for various 
initial conditions in figure 4. 
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5.5. Nil geometry 

When we take A = — 1, // = v = 0, G is isomorphic to the Heisenberg group. The flow 
equations 

dA . dB dC „ A n , . 

— = -8A 4 , — = AB A 3 , — = AC A 3 , 36 

dt dt dt ( J 

can be easily integrated to yield 

A= (24t + Ao 3 )- 1/3 , B = fiov 7 ^^ 172 , C = A- l ^/(B y/A^). (37) 

Therefore, A — > whereas .B and C diverge. In the terminology of [20], this degeneracy 
is called a "pancake degeneracy" . 

5.6. R 3 

In this case, A = [A — v — 0, the metric is flat and the flow equations are trivial. 

5.7. H 3 , S 2 x R, andH 2 x R 

These cases do not arise from the construction that uses a 3 dimensional Lie group as 
described, however the metrics can be written explicitly for each case. The metrics 
for the three geometries are of the form g = A(t)g H 3, g = D(t)g^ + E{t)gs* and 
g = D(t)crsi. + E(t)g H 2, where g^, gs 2 , 9r and gu 2 are the standard metrics on 
the hyperbolic 3-space, the unit 2-sphere, the real line and the hyperbolic 2-space, 
respectively. Each of the resulting metrics is already conformally flat, therefore they are 
fixed points of the Cotton flow. Note, however that the Ricci flow degenerates if 2 x M 
and S 2 xR [20]. 



6. Conclusions 

We have introduced a new flow in three dimensions, whose fixed points are conformally 
flat metrics. We have presented a gradient formulation of this flow by finding an entropy 
functional. We have also shown that the commutator of this flow with the Yamabe flow, 
which preserves the conformal class of a metric, is another flow preserving conformal 
class. We have studied the evolution of the nine homogenous geometries in detail. Four 
of these, R 3 , H 3 , H 2 x R, S 2 x R, are fixed under the flow. Note that the last two of 
these are degenerated by the Ricci flow. Every homogenous metric in the SU(2) class 
converges to the round S 3 metric. For the Isom(R 2 ), Solv and Nil cases, the equations 
can be solved analytically. Every metric in the Isom(R 2 ) class evolves to a flat metric. 
The metrics in the Nil class tend to a pancake degeneracy. The metrics in the SL(2, R) 
and Solv classes develop cigar degeneracies. 

There are several points which require further study. As in the Ricci flow case, one 
can define solitons: 

Cij + V,<, + V,0 + Xgn = 0. (38) 
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A metric satisfying this equality, which we call a "Cotton soliton", can be thought of 
as evolving under the Cotton flow by just a diffeomorphism and/or a scaling of the 
underlying manifold. In the case of a compact manifold, the conservation of total 
volume implies that A vanishes; however, to allow for non-compact manifolds, we keep 
it. We have not been able to find a non-trivial solution of this equation, or of the 
simpler gradient Cotton soliton equation where £ 3 - = Vj<£>. However, it may be useful 
to consider the techniques introduced in [251 126] in this regard. Another open question 
is to determine the conditions on the initial metric for which the short time existence 
problem can be solved, but this seems to be much more difficult than in the second 
order case. One obviously interesting avenue, for geometric purposes, is the study of 
singularities formed under the Cotton flow, or under other possible flows obtained by 
combining the Ricci and Cotton flows. It would also be interesting to find a nonlinear 
sigma model whose beta function is the Cotton tensor. 

Appendix A. Computation of the flow equations 

Here we show the details of the calculation of the Cotton tensor and how (|15p is derived 
using the differential forms. The connection 1-forms uo a & can be computed using the 
vanishing of torsion: de a + u a & A e b = 0. For the metric (fl3|) they specifically take the 
following form: 

u 1 2 = - 1 (XA + fiB- uC)e 3 , u 1 3 = - 1 (-XA + nB- vC)e 2 , 

u 2 3 = -^ = (-XA + uB + uC)e 1 . (A.l) 

Vabc 

Using the formula R a & = doj a & + oj a c A u c &, one computes the curvature 2-forms as 
R 1 2 = ~Xbc^ 2a2 + ~ 3u2 ° 2 ~ 2X V AB + 2XuAC + 2fiuBC)e 1 A e 2 , 

R 1 3 = —1—(\ 2 A 2 + v 2 C 2 - 3^ 2 B 2 - 2XvAC + 2X/J.AB + 2^uBC)e 1 A e 3 , 
ABC 

R 2 3 = J^ji^B 2 + v 2 ° 2 ~ 3X2 a2 ~ ^uBC + 2X^AB + 2XvAC)e 2 A e 3 . 

(A.2) 

Next, we compute the Ricci 1-forms and the curvature scalar using the formulas 
(Ric) a = LbR b a and 1Z = i a (Ric) a , respectively, as 

(ifcc)i = -^—(X 2 A 2 - fi 2 B 2 - v 2 C 2 + 2/2uBC)e\ 
ABG 

(Ric) 2 = _|_(/i 2 5 2 - v 2 C 2 - X 2 A 2 + 2XuAC)e 2 , 
ABC 

(R ic ) 3 = -J—^c 2 - X 2 A 2 - fi 2 B 2 + 2AM#)e 3 , 
ABC 

K = J^j(- x2a2 - V 2 B 2 ~ u2c2 + 2A//AB + 2XvAC + 2\lvBC). (A.3) 
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The Cotton 2-form is given by C a = dY a + u a & A Y b , where Y a = (Ric) a — jR-e a , which 
leads to 

C 1 = -^L_ [-A a A 2 fcB + uC - 2XA) - (jiB + vC)(jmB - vCf]e 2 A e 3 , 
C 2 = (AB ^)3/ 2 l-V 2 B 2 (isC + XA- 2/jB) - (uC + XA)(uC - \A) 2 }e 3 A e\ 

^ 3 = (ABcy/i [ ~ u2c2{XA + fiB ~ 2vC) ~ {XA + ^ B){XA ~ /ijB)2]el A el 

(A.4) 

Finally, using <9 t e a = *C a , one obtains (flol) . 

Appendix B. Coordinate invariant computation of the Cotton tensor 

It is sometimes useful to consider the coordinate free version of the computation above. 
For this purpose, we follow [27] and [22]. Define the (0, 2) tensor s by the formula: 

s = Ric-jg. (B.l) 

A Cotton 3-form on a generic D-dimensional manifold A4 can be defined by 

C(X,Y)(Z) = (V x s)(Y,Z) - (Vys)(X,Z), (B.2) 

where X, Y, Z are vector fields on Ai . Let us now return to the case of a homogeneous 
3-dimensional manifold. Recall {F±, F 2 , F 3 } from section [5j First, we compute s using 
the formulas for Ric(Fi, Fj) and TZ directly from [22J: 

S (F 1 , Fx) = -^—((X 2 A 2 - (fiB - vC) 2 ) - \(2/iuBC + 2vXAC + 2XfiAB 
BG 4 

- \ 2 A 2 - fi 2 B 2 - v 2 C 2 )), 

S (F 2 , F a ) = J^d^B 2 ~ ( XA - vC f) ~ \{^vBC + 2v\AC + 2XfiAB 
-X 2 A 2 - l i 2 B 2 -v 2 C 2 )), 

s(F 3 , F 3 ) = ^—{{n 2 C 2 - (XA - fiB f) - -(2nvBC + 2v\AC + 2\\iAB 
AB 4 

- \ 2 A 2 - fi 2 B 2 - v 2 C 2 )) (B.3) 
and clearly, s(Fi, Fj) — if i ^ j. One has 

(Vxs)(Yi Z) = Z)) - s(VxV, Z) - s(Y, V X Z). (B.4) 

When {X,Y,Z} C {F!,F 2 ,F 3 }, the first term on the right is by homogeneity. 
V ' FiFi = for any i, and 

-\A + jiB + vC 



V Fl F 2 = -\[r 1. r 2 \ - \<i(ir L ) r 2 - \ <i<t r 2 1 r X ] = ^ ' 

_ ,. f , -i/C + AA + y g 

Vf 2 ^3 — ^ ^l) Vfs-ti — — -r 2 , 





,F 2 ]-(adF 1 ) : 




+ i/C + XA _ 




A Fl 


//S - 


XA-uC _ 




C 


AA- 




B 



, vC-fiB-XA 

vf 2 F\ — -p, F 3 , Vf 3 F 2 — — Fi, 

V Fl F 3 = ~^ F 2 . (B.5) 
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From the definition, C(X, Y)(Z) is antisymmetric in (X, Y), we have C{F i) Fj)(Fj) = 
for every Since the other antisymmetry is not apparent, we prove: 

Lemma. For any C(Fi, Fj)(Fj) = 0. 

Proof: We may assume i ^ j since the equality case is taken care of above. Let k be 
the third index, different from i and j. Recall that s is diagonal. Now, 

tf(F 4 , = (V^Xi^F,) - (V^XF,,^) 

= s(V Fl F 3 , Fj) - s(Fj, V Fl F 3 ) + s(V Fj Fi, Fj) + aft, V F] F 3 ) 
= s{F j ,aF k )=0, (B.6) 

where a is some constant. Note that we have again used homogeneity to drop terms 
that involve d F .. This finishes the proof. 

Therefore, the only interesting components of the tensor are the ones involving all 
3 indices: 

C(F U F 2 )(F 3 ) = -s(V Fl F 2 , F 3 ) - a(F 2 , V Fl F 3 ) + s(V F2 F u F 3 ) + sft, V F2 F 3 ) 

= s C XA + » B + UC F s , F 3 ) - a(F 2 , AA -^-^ F 2 ) 
+s{ ,B-XA-vC ^ + s(Fi) -,B + ,C + XA Fi) 

\ XA — vC — uB . ^ . -aB + vC + XA . ^ 

= -2us(F 3 , F 3 ) ^s(F 2 , F 2 ) + aft, Fj, 

and the other Cfti, Fj)(F k ) can be computed in a similar way. Defining C^k = 
C(Fi, Fj)(Fk), the flow equations ( |T5l) can be written as: 

In fact, in a general coordinate frame, to relate this formulation to the one presented 
in the main text, we define 

Qjk = C(di,dj)(d k ) 

= (V l s)(d 3 ,d k )-(V 3 s)(d l ,d k ), 
= di(s(dj, d k )) - s(Vidj, d k ) - s{d jl VA) 
- dj(s(di, d k )) + s(Vjdi, dk) + s(di, W 3 d k ) . 
Using s(Vidj, d k ) = s(T e i3 dg, d k ) = Y l y s(di, d k ), this simplifies to 

Cij k VjSjfc Vj'Sjfc, 

where Sj k = s(d 3 ,d k ). In the case of 3 dimensions, one can get a symmetric, traceless 
and covariantly conserved 2-tensor out of this by contracting it with the completely 
antisymmetric Levi-Civita tensor: 

C 1 — ^ knm s~i 

2y/g 
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